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Abstract 



In a nonstandard way we split up the polarized quark distributions qt and 
Hl into their twist-two, quark-mass, and interaction-dependent parts, empha- 
sizing the sensitivity to quark intrinsic transverse momentum. We show how 
to derive the Burkhardt-Cottingham sum rule in this approach and derive 
a similar sum rule for the chiral-odd distribution /12. The effect of intrinsic 
transverse momentum in experimental observables is illustrated in the calcu- 
lation of the 0(1/(5) double-spin asymmetry Alt in Drell-Yan scattering. 
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I. INTRODUCTION 



Stimulated by recent |I| and planned ||^ experiments with polarized electron and proton 
beams, spin-dependent distribution functions have received a lot of attention. Due to increas- 
ing accuracy, measurements of higher-twist distributions seem feasible in the near future. 
They are particularly interesting, since they contain valuable information on quark-gluon 
correlations. The goal of this paper is to analyze the role of quark transverse momentum 
in the twist-three distributions^ gri^) and or equivalently the linear combinations 

g2{x) = grix) - gi{x) and h2{x) = 2/ii(x) - 2hi{x). 

At leading-order level, one may consider transverse-momentum-dependent distributions 
which can be measured in Drell-Yan (DY) scattering at small but fixed Qt In deep 

inelastic lepton-hadron scattering (DIS) or in DY integrated over Qt one measures the 
distributions integrated over transverse momentum; these are the twist-two unpolarized 
momentum distribution fi{x) and two twist-two spin-dependent distributions, the helicity 
distribution gi{x) and the transverse spin distribution hi{x). At leading order it does not 
matter what the transverse momentum dependence is; the observables (i.e., structure func- 
tions and asymmetries) do not change if one assumes zero transverse momentum [i.e., 
distributions. At subleading 0{1/Q), the situation is less simple. As was demonstrated in 
Ref. [Q, for instance, assuming free quarks without transverse momentum gives the non- 
physical result g2{x) = 0. Thus, even though g^ and do not depend explicitly on k^, 
assuming zero transverse momentum may considerably alter their shape (in x). We will 
make these statements more precise in due course. 

First, let us recall some facts about the transverse spin-dependent distribution grix) 
I^J^. Using operator-product-expansion (OPE) techniques, one shows that grix) can 
be decomposed in three pieces, a twist-two part depending on gi{x), also known as the 
Wandzura-Wilczek (WW) term ||^, a quark-mass part depending on hi{x) P, and an ex- 
plicitly interaction-dependent part. The latter is the most interesting, since it is sensitive 
to off-shellness and confinement effects. Instead of using the OPE, we rederive the different 
terms, using nonlocal matrix elements only. In this approach the transverse momentum of 
quarks is an essential ingredient. The nonlocal matrix element that fixes the interaction- 
dependent part contains two good quark fields and one transverse gluon field, so that it 
can be given a clear parton interpretation 0. Of course, if moments are taken, one should 
recover the OPE local matrix elements. In this formalism it is also straightforward to derive 
the Burkhardt-Cottingham (BC) sum rule Jq dx g2{x) = [0], although one must make 
assumptions. The status of the BC sum rule is less rigorous than, for instance, the Bjorken 
sum rule. It was claimed that in perturbation theory it is valid at least up to 0{as) ||11|| , 
although recently this claim was disputed 



The twist-three distribution function hiix) [13] is less well-known. It does not occur in 



To avoid confusion, we emphasize that we are sometimes using names reserved for the structure 
functions measured in inclusive deep-inelastic lepton-hadron scattering for quark distributions. 
In fact, the distributions should be carrying a quark (or antiquark) flavor index, which is often 
suppressed to simplify the formulas. At the end the structure functions are obtained as a simple 
weighted sum over quark distributions. 
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the inclusive DIS cross section other than multiphed with a quark-mass factor (hke hi{x), it 
is chirally odd). Nonetheless, one can analyze it using the OPE-based techniques [|14|. Again 
one finds three parts; a twist-two part depending on hi{x), a quark- mass part depending 
on gi{x), and an interaction-dependent part. The interaction dependent part can also be 
given a parton interpretation [|1^]. Using the same techniques as for qt, we derive within 
our approach these parts, expressing them in a finite number of nonlocal matrix elements. 
We derive a sum rule very similar to the BC sum rule, namely, Jq dxh2{x) = 0. 

In order to support our claim that explicit treatment of quark transverse momenta is 
important at twist-three, we consider Q^-averaged polarized DY scattering. Jaffe and Ji 
derived an 0{1/Q) asymmetry, Alt, for longitudinally-transversely polarized DY scatter- 
ing [Q. In their calculation they assumed zero quark transverse momentum. Without this 
assumption we find modifications in the asymmetry. To estimate the difference between the 
two results we use bag-model distributions. The difference turns out to be considerable. 

This paper has the following structure. In Sec. y, we analyze the content of grix) and 
using Lorentz symmetry, discrete symmetries and the QCD equations of motion. We 
decompose the distributions in the three different parts, and derive the sum rules. In Sec. |T| 
we calculate the polarized DIS hadron tensor, emphasising the role of transverse momen- 
tum. In Sec. |V|, the more involved but essentially similar hadronic tensor for polarized DY 
scattering is considered, from which follows the double-spin asymmetry. 



II. ANALYSIS 

In this section we decompose the polarized twist-three distributions grix) and hi^x) 
into their building blocks of different dynamical origin. Instead of using conventional OPE 
techniques, we work directly with the nonlocal quark-quark and quark-gluon-quark matrix 
elements. In the calculation of a specific cross section these matrix elements appear as the 
non- calculable parts in the diagrammatic expansion of the amplitude, i.e., those parts that 
connect hadron lines to quark or gluon lines. 



A. Correlation functions 

First, we will expand the quark-quark and quark-gluon-quark correlation functions as 
needed up to 0{1/Q), using hermiticity, as well as Lorentz invariance, parity invariance, 
and time-reversal invariance. 

The correlation function describing a quark of fiavor a in a spin- 1/2 hadron A, is defined 

by m 

i<^a/AUiPS;k) = I ^,e^'-''{PS\^^;\0)^pi^\x)\PS), (2.1) 

diagrammatically represented by the blob in Fig. |]. Here P is the hadron momentum vector 
(P^ = M^), and S is its spin vector {P ■ S = 0, and S"^ = —1). In lightcone variables (for an 
arbitrary four- vector p = [p~ ,P'^,Pt], with = (p*^ ±p^)/\/2) they read 
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p 

s 



2P+ 



-A- 



M 
'2P+ 



P+ 



(2.2) 
(2.3) 



where A is the hehcity, and St = {S^, S"^) is the transverse polarization. Note: A^ + S"^ 
— S"^ = 1. The quark momentum is written as 



2xP+ '■ 



(2.4) 



with the hghtcone momentum fraction x = k^/P^. 

An important step in our analysis is the expansion of the correlation function in terms 
of all possible Dirac matrices multiplied by scalar functions depending on k ■ P and k"^, and 
vanishing when these variables become larger than a characteristic hadronic scale A^. Using 
the demands of hermiticity, Lorentz, parity, and time-reversal invariance, one finds |^,^ 

^{PS; k) = ail + (as/M) f + (ag/M) 

+a4 75)^ + (a^/M) 75 [f, ^] + {a,/M) -f,[^, ^] 

Hay/M') {k ■ 5)75^ + (ag/M^) {k ■ 5)75^ + (ag/M^) {k ■ S^^, (2.5) 
where the dimensionless amplitudes aj(o", r) are real functions of the scalar variables 

2k ■ P 2k- P^ 



a 



T = 



M2 



X- 



M2 
2k-p+ 
M2 



+ X, 



(2.6) 
(2.7) 



Since the relevant object in deep inelastic processes is the integral over k of $(PS'; k), it is 
convenient to define the projections (suppressing flavor and hadron labels) 



dk-i:i[T^{k)] 



1 f dx d Xt 
2ti (27r)2 



exp[i(xP+x- - kr ■ xt)] {PS\i}){{))Ti){{),x- ,Xt)\PS), (2.8) 



with r taken from the complete set {1, 275, 7^", 7^75, ia^'^'j^}. Consider for instance the 
projection for F = 7"'": 



ni^] = ^Jdadr6{ 



fc2 

r - X(T + + j^) [02(0-, r) + xa3{a, r) 



(2.9) 



where we inserted the expansion of the quark correlation function in terms of amplitudes, 
Eq. (|2.5|) . Clearly, it is an 0{1) function, /i, depending on x and fc-r- Similarly, for F = 1 
one has 



m 



M 1 

P+M2 



r fc2 

/ dadr 6{t — xa + x'^ + -Jr-) ai (a, r 



(2.10) 
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which is again an 0{1) function, e, of x and k'^, but multiphed by a factor M/P~^. In deep 
inelastic processes, characterized by a large momentum scale Q, this factor will give rise 
to an M/Q suppression. The following projections are leading (here and in the rest of the 
paper, the indices i and j are transverse, i.e., taking the values 1, 2): 



$[7+] = /i(x, 4), 

"^"[7^75] = 9il{x, k%)\ + ^it(x, fc|) ^ , 
^[iiy'^^lb] = hiT{x, k'^)S'rp + hj^^ix, k^)X + h^'^i^, k%' 

where we have defined 

If k"^ 

k%) = — j dadT5{T - xa + x"^ + j^) [02 + xa^] 



k 



M 



r\jrp 



(2.11a) 
(2.11b) 

(2.11c) 



giLix, kl) 
girix, kl) 
hirix, fc|) 



M2 

1 

Jp 
1 

AP 
1 
M2 

1 

AP 
1 
M2 



fc2 

dadr 5{t — xa + x^ + j^) '^4 + (^o" — a;) (07 + xa^) 

fc2 

dodr 5{t — xo + x"^ + j^) [—<^7 — xag] , 

i f k^ 

— J dadrdir - xa + x"^ + j^) [-2a5 - 2x06] , 

fc2 

dadr 5{t — xcr + + -^^) 2a6 — 2(|o" — a;)a9 

fc2 

dadr 5{t — xcr + + -^^) pag] . 



(2.12a) 
(2.12b) 
(2.12c) 
(2.12d) 
(2.12e) 
(2.12f) 



These leading transverse momentum distributions can be directly observed in DY at mea- 
sured Qt [0. For DIS or Qr-averaged DY one needs only specific fc-p- integrals of Eqs. ( ^.lip . 
Particularly, the ones occurring at leading order inl/Q are 



j rf2fcT$[7+] = / d^krhix^kl) ^ /i(x), 

j rf2fcr$[7+75] = d^krgiLix^kl) = \gi{x) 



c/2fcy$[i(T*+75] = Sir / d^kr 



hirix, k^) + 



2M2 



= Slhi(x) 



(2.13a) 
(2.13b) 

(2.13c) 



Note that gir and hj-^ vanish because their fcr-structure is odd (under the operation kx 
— fc^^). We recognize the longitudinal momentum, helicity, and transverse spin distributions. 
Expressed in terms of the amplitudes, with the help of Eqs. ( |2.12| ), they are 



/i(x) = dadT9(xa — t — x^) [02 + xa^] , 
gi{x) = dadrO^xa — t — x^) 04 + (|cr — x)(a7 + xag) 
hi(x) = n dadrO^xa — r — x^) — 2a5 — 2xa6 + {xa — t — x^)a9 



(2.14a) 
(2.14b) 
(2.14c) 
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The ^-function together with the spectral condition {P — k)"^ > leads to an integration 
area as illustrated in Fig. Note that for x = 1 the integration area becomes zero and 
consequently the distributions vanish. 

For the subleading projections, we will need only the integrals over kx- Three Dirac 



projections remain |n 



w 

p+ 

~M 
P+ 



J d'^kr^l] = e{x), 
Jd'kTnYi,] = Sir9Tix), 
^ j d^kT<^[ia+-^,] = \hL{x). 
These twist-three x-dependent distributions in terms of the amplitudes read 

e{x) = dadrO^xa — r — x^) [ai] , 



grix) = IT dadrO^xa — r — x^) 



0-4 — |(xcr — T — x^) 



hL{x) 



TT 



dadr 6(xa — t — x 



-2a5 — eras + 2(^0" — x^ai^ 



(2.15a) 
(2.15b) 
(2.15c) 

(2.16a) 
(2.16b) 
(2.16c) 



The chirally odd distribution e(x) contributes at O^i/Q"^) in unpolarized Drell-Yan scat- 
tering |1^. The polarized distributions grix) and hi{x) will be analyzed further in the 
next subsections. In summary, the quark correlation function where the nonlocality is pure 
lightlike (and in the — direction) can be parametrized up to and including twist three as 

2 [ d^krdk-^ik) = fi{x) 7~ + \gi{x) 757" + 5^/ii(x) i-f^a- 



+ 



M 
P+ 



e(x) 1 + S^^grix) j^-fi + Xhiix) i-f^a^ 



+ higher twist. 



(2.17) 



At 0{1/Q), another type of correlation function needs to be considered, namely, the 
quark-gluon-quark matrix element 



{Ma/A)ip{PS;k,p) 



d'^x d^x' 
(27r)4 (27r)4 



^ik-x i{p—k)-x' 



{PS\^l^y{Q)gAUx')^Pl-\x)\PS), (2.18 



pictorially represented in Fig. The demands of hermiticity, parity, and time-reversal 
invariance lead to the constraints 



[M'{PS;k,p)f = -f°M'{PS;_p,k)-f^ 
M'{PS;k,p) =-f'^Mi{P-S;k,P)l'' 
[M'iPS;k,p)]* = j,CMi{PS; k,P)Ch^o 



[Hermiticity] 

[Parity] 

[Time reversal] 



(2.19) 



where C = ^7^7° is the charge conjugation matrix, and k^ = k^. In deep inelastic processes, 
at least on the tree level, one is sensitive to the projections 
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M'[r]{x, kT) = \j dk- d^pTi \rM\k,p) 



= \j exp[z(a:P+x- - kr ■ Xr)] (P5|^(0M^(0)rV^(0, x", a;r)|P5), (2.20) 

specifically, at the 0{1/Q) accuracy level, to the leading ones with V = 7+, 7^75, or Z(T*"'"75. 
If Qt is not measured, then at 0{1/Q) one may also integrate over kr- Using the parity 
and time-reversal constraints of Eq. (|2.19|) , we define the following real quark-gluon-quark 
distributions: 

S't~9t{x)^ 2^ / d'^TMAgh^l, + i4l^] + h.c. 

= 4^ / {PS\m9ATm[9h''l^ + ^4l^m^^)\PS) + h.c, (2.21) 

XhL{x)= j d^kTM'[gTijia^^-i5] + \i.c. 

= 4^/ {PS\mgMmT^,^a^^lm^-)\PS) + h.c., (2.22) 
where iIj{x~) denotes ip{Q,x^ ,Qt) and we have used the tensors 

gl^" = g^"' - n^ln"^ - n\nt, (2.23) 
= e^%,<n-. (2.24) 

The lightlike vector n+ and n_ have only a non-zero + and — component, respectively, such 
that n+-n_ = 1. The reason for the suggestive nomenclature of grix) and hL{x) will become 
clear later. 



B. Decomposition of grix) 



The first important step in the splitting up of gxix) into its constituent parts is the 
implementation of the QCD equations of motion for the quark fields ipip = mip In 
lightcone language this implies the elimination of the 'bad' field = A^ip = (7~''7~/2)-?/' 
occurring in grix) [Eq. ( p.l5b|) ] in favor of the 'good' fields 4'+ = ^+4' = and 
A^. Explicitly one has 







d'^x 
(2^ 
d^x 

d'^x 
(2^ 



gifc.x^p^l^^Q) i/)(0) -m ia'+-f5ij{x)\PS) 



Ak-x 



{PS\ij{<d)\]/i + g - m]ia'+-i^^{x)\PS) 



e^^-^(P5|^(0) [fc+7- + (4 + g^Ti!^))l3 - ia'+^^'4){x)\PS) 



(2.25) 
(2.26) 
(2.27) 



where translational invariance of the hadron states is used and the lightcone gauge = 
is chosen. Performing the integral | / dk~ , one obtains 

= A;^$[7+75] + ^e^kT,<^>b+] + m$K+75] + M,[gi^j+j, + ze^7l, (2-28) 
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using the definitions ( p.8|) and ( |2.2CI| ). Finally, upon taking the real part0, integrating over 
and using the parametrizations for the projections given in Eqs. (|2.11bD , ( p.l3c|) , (|2.15b ), 
and ( |2.21| ), one gets 



9t{x) 



2M2 



+ 



X 



M 



X 



(2.29) 



This result will be used in the next sections, when we will calculate DIS and DY structure 
functions. We make a few remarks. First, at 0{1/Q) the transverse covariant derivative 
inevitably introduces a sensitivity to the transverse-momentum-dependence comprised in the 
first term on the lefthandside of Eq. (|2.29| ), even if transverse momentum is not observed 
directly. Secondly, this term contains the functions qit, to which one is not sensitive at 
leading order if transverse momentum is not observed [see Eq. ( p. 13b )]. If this distribution 
is proportional to a delta function 5(fcy), the integral vanishes. So if one assumes zero 
transverse momentum, the difference between grix) and grix) is of order m/M. In general, 
however, this cannot be assumed. 

It is possible to re-express the above relation in terms of the fc^^-integrated distributions 
only. For this one needs to use the expressions in terms of the amplitudes aj(cT, r) given in 
the preceding subsection and the following relation for a linear combination of amplitudes 
J^{x,a, t): 



— TT j dy j dadr 6{ya — T — y'^ 



[a - 2y)J^{y, a, r) + (ya - t 



y )"^(2/,c^,i") 



(2.30) 



which is most easily proven by differentiating both sides with respect to x and using the fact 
that the integration area vanishes at x = 1. Applying it to the term containing girix, k"^) 
in Eq. ( |2.29| ), using its expansion in amplitudes Eq. ( p.l2c| ), this gives after some algebra 



9t[x) 



dy h — 



y 



M 



hi{x) 



X 



'dy'^^y^ 



+ 9t{x) - / dy , 

y 



(2.31) 



The first term was derived by Wandzura and Wilczek in The quark- mass terms are 
chirally even, since both m and hi{x) are chirally odd. The last two explicitly interaction- 
dependent terms only need the nonlocal matrix element grix), defined in Eq. ( |2.21| ). 
With the relation ( |2.31| ), one easily checks that 



dx g2{x) 



dx [grix) -fi'i(x)] 



r dx G{x) - r dx —G{y), (2.32) 

JO JO Jx V 



with G{x) = —gi{x) + {m/M)[hi{x)/x] + grix). Provided that in Eq. ( p.32| ) one may 
interchange the order of integrations, this amounts to the Burkhardt-Cottingham sum rule 
jQdxg2{x) = ig. For more details concerning its validity, see Refs. P,[5|,P^ . 



^Note that ^[T] is real if F is taken from the set {1, i^^, 7^, 7'^75, ia'^'^^^}. 
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In Ref. Jaffe and Ji calculated the contributions of different dynamical origin in grix) 
for a massless quark in a nucleon bag. In Fig. H their results are plotted. Also plotted 
is gri^x) [the authors call this function U{x)]. Note that, as we argued before, the role of 
transverse momentum is fairly important, because else gr^x) would equal grix), as follows 
from Eq. ( p.29[ ) with a 6{k'^) dependence and m = 0. Nonzero transverse momentum is an 
important feature for a confined quark. Note finally that the BC sum rule in the bag reads 
J'^^dxg2{x) = 0, because violation of translational invariance leads to a support between 
— oo and oo. See also Fig. |. 



C. Decomposition of hiix) 



The twist-three distribution hiix) can be treated in a similar fashion. This time, 
implementing the QCD equations of motion, one must consider the matrix element 



(P5|V(0)[ip(0) - m\j^j5^{x)\PS) -- 
This leads to the relation 

m gi{x) 



hiix) 



M2 



X 



M X 



0. 



+ hLix) 



(2.33) 



(2.34) 



involving the transverse momentum distribution h^^{x,k'^) to which one is not sensitive 
at leading order [Eq. (|2.13c|) ]. The fc^-integral can be re-expressed as an x-integral by 
means of the relation ( |2.30 ) using the amplitude expansion Eq. ( |2.12e|) . We obtain for the 
decomposition of hiix), 



hhi^x) = 2x 



hi{y) 

dy — — 

y 



m 
M 



9i{x) 



2x 



X 



dy — — 

yi 



+ h^ix) — 2x1 dy 



hiy) 



(2.35) 



The compact expression of the interaction-dependent part in terms of /^^(x), defined in 
Eq. (|2.22|) , may be a practical one for investigating models. Relation ( p.35| ) was also derived 
by Jaffe and JiQ using OPE techniques |p!^ |. 
With the above relation in hand, we consider 

[ dxh2{x) = 2 I dx [/iz,(x) — hi{x)] = [ dxH(x) — [ dx2x [ -^H(y), (2.36) 
Jo Jo Jo Jo Jx y^ 

with H{x)/2 = —hi{x) + {m/ M)[gi{x) / x\ + hiix). Provided that the order of x- and y- 
integration of the second term may be interchanged, we find the sum rule Jq dxh2{x) = 0. 



^ We find a discrepancy of a factor of 2 in the third terms on the righthandsides of Eq. ( p. 35 ) and 
Eq. (56) of Ref. [0] . The lower integration limit and the power in the integrand in the latter case 
are obviously misprints, Cf. Eq. (51) of the same paper. During their calculation, Jaffe and Ji also 
introduce a function hi[x), which is not the same as ours. Because in their case it is an auxiliary 
function, which doesn't occur anymore later, and as good names are scarce, we considered it safe 
to re-use the name. 
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The validity of this sum rule, of which no reference is known to us, crucially depends on the 
small-x behavior of the function H{x). 

In Ref. the different contributions to /^^(x) were calculated for a massless quark in the 
bag. Their results are plotted in Fig. A comparison with the bag- model /^^(x) in the same 
figure shows that the zero-quark-mass zero-transverse- momentum relation /il(x) = hiix), 
following from Eq. ( |2.34|) , is severely violated. Remarkably, h2{x) = 2g2{x) in the bag (see 
Fig. ^), and it satisfies the above sum rule, provided that one extends the integration region 
to the full line f— oo, oo). 



III. POLARIZED DIS 



In the following two sections we address the question how the different distributions 
of the preceding section occur in the physical observables. In this section the polarized 
deep inelastic cross section through order 1/Q pO| is calculated as an illustration of the 
more complicated DY calculation. The diagrammatical method by Ellis, Furmahski, and 
Petronzio |^ is used. 

First, let us discuss the kinematics. The target hadron has momentum and spin vectors, 
given in lightcone representation in Eqs. ( p.2|) and (|2.3| ). The incoming virtual photon is 
assumed to have momentum q such that —q^ = is large, but Xbj = Q^/2P ■ q is constant, 
and can be chosen^ 

4- 1 



2xbjP+ 



(3.1) 



Instead of working with explicit vectors P and g, it is often more convenient to express them 
in terms of the lightlike vectors n+ and n_, defined as 



V2 

[0, /«, Ot] ~ -Q^bjP, 



K 



V2 

— {q + XbjP), 



(3.2a) 
(3.2b) 



where the parameter k = \^XbjP^/Q- 

Consider the quark Born diagram of Fig. |^ (the antiquark diagrams can be obtained by 
crossing). It stands for 



(3.3) 



Where the sum runs over quark flavors a with charge in units of e. Using the lightcone 
representations for k and q, it is easily checked that 



SUk + qf 



m 



5{k+ + q+) 
2q- 



(3.4) 



^Here and in the following we will denote by '~' that we have neglected 0{1/Q'^) contributions. 
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equating x to Xbj- Furthermore one can write 

/(f + ^ + m + fcr7j + (3-5) 

Note that the only ^"-dependence of the integrand in Eq. ( pl3|) resides in $(/c), hence indeed 
we are sensitive to / dk~ ^{k). From exphcit calculation (or from dimensional arguments) 
the result is obtained from the first term in Eq. ( p. 51) resulting in an expression that 

only contains the leading projections $[r] of Eqs. ( |2.11| ) integrated over fc-r- Therefore, at 
leading order, one is sensitive to the (chirally even) distributions fi{x) and gi{x), defined in 
Eqs. (g^). In fact, 

2M W^^, =Y.el [-9't mx,,) + 9ti^t,)] , (3.6) 



where the transverse tensors are defined in Eqs. ( |2.23| ) and ( [2.241 ). This leading order result 
is electromagnetically current conserving, that is ~ 0- 

Turning to the 0{1/Q) contributions, we find that they either come from the first term 
in Eq. ( |3.5| ) combined with the subleading projections, or from the combination of the last 
two terms in Eq. ( p. 51) and the leading projections. For the former combinations the same 
argument holds concerning the fcj^- integral, hence we only need the x-dependent distributions 
of Eqs. (|2.15| ). For the latter combinations the argument goes for the m-terms, but the 
second term of Eq. ( ^.5] ) selects the transverse momentum distributions that multiply a 
kjn-odd structure, i.e., giT{x,k'^) and hj-j^{x, k"^). Thus, after performing the trace and the 
integrations, we end up with the (9(l/(5)-part of the Born diagram 

M 



2Miyr 



Bl 



-V2i 



f 2 a 2 

J ^ ^'^2J^^'^Ti^bj,krp) 



—h^(x,- 



where Sj^ 



9t 



(3.7) 

by means of the 



We obtain a more physical picture if we eliminate grix) 
relation (|2.29| ) found in the preceding section, so that the result contains only distributions 
which are matrix elements of the independent 'good' fields and A^- We get 

., , ... „ 



2M W^'" 



Bl 



Q 



Y^el{2te^%,S!^r[J d'k 



+ 



m 



M 



paSrp'n_ Xbjg'jri^xifj 



(3.8) 



where q = q/Q ^ (n_ —n^)/\/2. Note that the first two terms are current conserving in the 
sense that they give if contracted with q. This is as expected, since for the theory without 
interactions, i.e., without gluons, one has grix) = 0, and the Born diagram is the complete 
result. In QCD, however, one must include the gluon diagrams of Fig. ^ which contribute 
at 0{1/Q) as well. Written out in full, they are 



2MW, 



Y^elj d'kd'p5{{k + q) 



2 2 

m 



x<^Tr 



-Tr 



«7i 



+ ^ + m) 7 



+ d - m 



7 



(3.9) 
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The delta function can be approximated as in Eq. ( p.4| ). The extra parton momentum again 
has p^^Pt -C Q^. At 0{1/Q) accuracy, the polarization sum may be replaced by q~'~^^. 
Realizing that 7'''7i7''' = 0, one can effectively use for the hard propagator 



1 



7 



2g- QV2' 



(3.10) 



Since this does not depend on p anymore, the corresponding integral works on the M 
directly and we are sensitive to the projections (p.20|) or the complex conjugates of them 



(using the hermiticity condition ( |2.19| ) to exchange the order of the arguments). Particularly, 
in leading order, one has only the projections involving 7"*", 7^75, or ia'^^'j^. Also, in leading 
order, one only needs the fcr-integrated projections. Finally, it turns out that exactly the 
combination (|2.21|) , parametrized by grix), occurs; 



As expected, this combines with the Born diagram ( p.8|) , such that the total (9(l/Q)-result 

2M 



2M«r + W^J 



uu 



e4I 



2xbjM 



(3.12) 
(3.13) 



is current conserving. Observe that the final result can be expressed solely in terms of the 
distributions g^{x) = gi{x) + (^) turn out that the DY-case is not that simple. 

The flavor sums in the zeroth order result Eq. ( |3.6| ) and the first order result Eq. (p.l3|) 
may be extended to include antiquarks, provided that one replaces /° — *■ and g"" — > —g"", 
in accordance with charge conjugation invariance. Using the standard decomposition in 
structure functions one sees that Fi{xbj) = {l/2)J2a^a[fi{^bj) + fi{xbj)], and gi{xbj) = 
(1/2) Ea el [gtixbj) - gK^bj)] with z = 1, 2. 



IV. POLARIZED DY 
A. Kinematics 

The DY hadron tensor is a bit more involved, since it contains two hadronic blobs. 
The principle steps, however, resemble closely those of the preceding section. First, some 
kinematical preliminary remarks. In the Drell-Yan process or massive-dilepton production, 
A + B ^ i + i + X , the two spin-^ hadrons have momenta Pa and Pb, respectively, and 
are assumed to be pure spin states. That is, their spin vectors satisfy S\ = = — 1. 
The measured lepton pair with momenta ki and k2 has total momentum q = ki + k2- The 
relevant lightcone-coordinate representations are listed below. 
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Pa = 




Ot 


) 


Pb = 


\p- 


Ot 




OA — 






s 1 


Sb = 




. Mb ^ ' 


Q = 


XbPb,Xa 


P+ 





(4.1) 

(4.2) 

(4.3) 

(4.4) 
(4.5) 

We work in the Drell-Yan limit where = and s = {Pa + Pb)^ become large with fixed 
ratio r = Q'^/s. Also we consider only dileptons with transverse momentum = ^ A^. 
We again define lightlike vectors n+ and as in Eqs. (^, with k = v^x^P^/Q ~ 
Q/{^/2xbPb)- The transverse tensors (/j.*^ and et^ are defined in terms of them as before 
[Eqs. ( p.23| ) and ( p.24| )]. Important is to note the precise meaning of transverse here. A 
transverse vector aj. = gl^'^ai, has both ■ Pa = and aj- ■ Pb = 0. In general the photon 
does have transverse components. 

The dominant elementary process underlying the reaction is the annihilation of a quark 
(antiquark) of hadron A by an antiquark (quark) of hadron B into a massive photon of 
momentum q which subsequently decays into the lepton pair. The cross section can be 
written as 



da 



a 



d^qdVt 2s Q 



-L W^'' 



(4.6) 



where the hadron tensor 

W'^-'iq^PASAlPBSB) 



d'^x 



(2vr)4 

is contracted with the lepton tensor 

L'^^ = 2 /tf A;^ + 2 ki^k'i - 



e''-^{PASA\ PbSb\ [^^(0), r{x)] \PaSa\ PbSb)- (4.7) 



(4.8) 



Since the latter is symmetric, we will henceforth discard the antisymmetric part of the DY 
hadron tensor. One can rewrite the lepton tensor as L^^ = —Q'^{g^'^ — q^q^ + IH^), where 
I = [ki — k'^IQ, satisfying / ■ g = and r = —1, defines the lepton axis. The angles d and 
are those of the lepton axis measured in a particular dimuon rest frame O' defined by 
Collins and Soper p2|. The polar angle is defined with respect to 



Pu-Pa ^ Pa- P. 



B 



^Bi 



(4.9) 



or rather its normalized version = j \l —Z'^ . Note that g ■ Z = and that Z does not 
have transverse components. The azimuthal angle is fixed with respect to q^ = qr — istr ■ 
with q^ = q^/Q. The four- vector q± is an example of what we call a perpendicular vector [Q. 
Defining 
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(4.10) 



a perpendicular four- vector a^^ = g^axv satisfies a_L ■ g = a_L ■ i = 0. For completeness we 
introduce 



The difference between 'transverse' and 'perpendicular' is of 0{1/Q), 

9t ^ 9± H 



Q ' 



n+ 



71- 



1 



q± 



ay ^ a_L H ^ — g, 

qt ■ bx ^ a± ■ b±, 
provided that ot • gr, ■ qr = 



(4.11) 

(4.12) 
(4.13) 
(4.14) 

(4.15) 

(4.16) 
(4.17) 



B. Hadron tensor 

We first calculate the DY quark Born diagram of Fig. ^. It gives 

Wb=\T. ^b-acl J d^ka d% 5\ka + h - q) Tr [^^{ka) ^ M^) Y] , (4.18) 

a,b 

where the antiquark correlation function $a is defined as in (O), but with ^ ij^^K The 
factor 1/3 is a color factor. The quark and antiquark momentum vectors read in lightcone 
coordinates 



kn 



k^ 4- h"^ 



(4.19) 
(4.20) 



and we assume that fc^, fc^j., fc^, and fc^^ all are of 0{\). Since then k^ <^ and k^ <^ , 
one can approximate the delta function, 



5\ka + h-q)^ 5{kt - g+)5(fc,- - g-)52(fc„T + k,T - Qt), 



(4.21) 



and Eq. ( [4.1^ ) reduces to 
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~ ^ E ^b-ael J d^KT d'hr S\Kt + hr - Qt) 

a,b 

y dk-^M) r (/ dk+Mh)) 7 



xTr 



(4.22) 



rendering Xa = xa and Xb = xb- The leading order result comes from inserting projections, 
Eqs. (|2.11| ), and the corresponding ones for hadron B, which can be obtained from the 
former by replacing + ^ — , f"' ^ f"', g"' ^ ~g°') and h"" ^ h°' |Q. One finds 



= -\T. ^b-ael I d^Kr d'hr 6\Kt + ^^r - Qt) 

$^[7+] $6 [7-] + $a[7^75] Ml'l^] 



9t 



+<i>4za^+75] $4z(T-'-75] {grr^^gA - 9T^Jg^'') > + 



(4.23) 



This expression still has the full dependence on Q, xa, xb, and Qt- If one does not measure 
Qt, one has to take the gr-integral of the cross section ( |4.6| ). In the dimuon rest frame O', 
the lepton axis vector / has only spatial components, and does not depend on qt anymore. 
Neither does q, so the qy-integral may be pulled through the lepton tensor to work on W^" 
directly. We will call the result after integration W^'^ . In the frame O', the components of 
Pa and Pb (or alternatively n+ and n_) have a q^- dependence. However, the vector Z does 
not. 

Returning to Eq. ( 4.2ij| ), we may in leading order replace gj!^ by g^^ according to 
Eq. ( [4.12| ). Integrating over qt in the dimuon rest frame O', the g±^s may be pulled 
outside the integral, since they are built from q and z [Eq. ( [4.10|) ]. Now the only qt- 
dependence resides in the delta-function, which is subsequently cancelled. The transverse 
momentum integrals then work on their corresponding projections directly, hence we may 
use the parametrizations (p.l3|). The result is the zeroth-order DY hadron tensor P| 



(4.24) 



= -\j:^l{ifi(^A)f!{xB) - \A\Bg1ixA)gKxB)]gT 

+ h1{xA)hl{xB) - iSA± ■ SB±)gT]} 



where we used the symmetrization of indices, S^j_S'^j_ = S'^j_S'^j_ + S'^^S'^ 



B±- 



Since in this paper we are interested in subleading order, we must find the 0{1/Q) 
content of Eq. ( [4.23| ). The corrections that we neglected above by replacing g^f by g'^ are 



'^'^Bl-l 



^ E ^ba el j d^kaT d^kbT 5'^ikaT + kbT - Qt) 



X 



$a[7^] $47 ] + '^'a[7^75] $47 75] 



+$ J^a*+75] <^b[ta'-l6]q^^ {g^^^L, + - gT^,q'2)]■ (4-25) 



Integrating over qt in the frame O' , we may pull q outside the integral. Aside from the delta 
function, the integrand is then linear in qt, selecting kT-odd structures in the projections, 
as can be seen from the resulting expression 
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U2 



Ma 



d'^kn 



.^g^^^{xA, kl) 



h\{xB) 



+XBq^^S''2^h1{xA) 



d^kr-^h^^ 



Mf 



IL 



(4.26) 



In addition, Eq. (|4.22|) contains contributions which are convolutions of a leading and sub- 
leading projection. The corresponding Lorentz tensors do not contain explicit ^T-dependence 
in the photon rest frame, so that the delta function vanishes and the kaT- and fc;,T-integrals 
work on the corresponding projections directly. Using Eqs. ( p.l3| ) and ( p.l5| ), we get 



+MA\A{q - XAhlixAMixB) 

-MA\B{q - XAg''T{xA)gl{xB) 

+MB\B{q + ^)^'^5a± Ki^A) XBhlixB)}. 



(4.27) 



If we eliminate grix) and /il(x) using relations (|2.29|) and (|2.34j) , respectively, like in the 



preceding section, we obtain W^^ = W^bi-i + 

= ^T.el{-XAZ^'S^J^ gtixA) \j d'k^^g^XB.k 



+Xaz^^S- 



B± 



Ma 



IL l^A, ^t) 



+\bz^^S''2^ 



d^kr 



-Xbz^^S"}^ hIixA) 



^diA^A, kl) 



h\{xB) 



9i{xb) 



M 



IL K^B, ^t) 



B 



+2m 



'Xbz^'^S'^^ hl{xA)gl{xB) - Xaz^^S]^^ gl{xA)hl{xB) 

-MB\A{q + z)^''S''b\ gt{xA) XB~gU^B) 

+MA\A{q - ^)^'^5b± XAhl{^A)h\{xB) 
-MA\B{q - zf^S% XAgUxA)9KxB) 

+MB\B{q + 5)^^5i hl{xA) XBhUxB)]. (4.28) 



As expected, the g-terms that violate current conservation only come multiplied with the 
explicitly interaction-dependent matrix elements grix) and hL{x). In QCD, one must include 
the gluon diagrams of Fig. |10|. Consider first diagrams (a) and (b) where the transverse gluon 
is inserted in hadron A, 
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a,b 



x|Tr 

+Tr 



p — q — m 
MHp, ka) YMh)il^ 



~ q — m 



■7 



(4.29) 



The delta function can be approximated as in Eq. ( [4. 211 ). Since the leading projections 
$6[r] are multipied by an overall 7"'", and since 7+7*7+ = 0, the hard propagator effectively 
reduces to 



1 



(4.30) 



Therefore, the p-integral can be pulled through the trace to work directly on the quark-gluon- 
quark correlations. Using the hermiticity condition ( |2.19| ) to relate M.{p^ka) to M.{ka.,p), 
one is sensitive to the projections ( p.20|) with F = 7+, 757''', or zcr*+75, or their hermitian 
conjugates. Working out the details, one find that to 0{1/Q) accuracy the Lorentz structure 
of the expression ( [4.29| ) in the frame O' has no qrj'-dependence, allowing us to use fc-r-averaged 
correlation functions in Eqs ( |2.21 ) and (2.22). The result is 



W 



gluon-A 



^ E el{-MAXA{q + z)^'S''J^ XAhUxAmxB) 

+MA\B{q + z)^^S'2^ XA~gUxA)g!{xB)}. 



(4.31) 



In the same way we find for the diagrams |T3(c) and |T3(d) with a transverse gluon insertion 
in hadron B: 



gluon- iJ 



-MBXB{q-z)^^S-2^ h1{xA)xBhl{xB)}. 



(4.32) 



The antiquark diagrams can be obtained from the corresponding quark diagrams by 
replacing m — *■ —m, e° — *■ —e"', — >■ g"" — and —>■ h"". However, the quark 
results are invariant under this operation, so one can use them, extending the flavor sum to 
include antiquarks. 

The total 0{1/Q) polarized DY hadron tensor contains no current non-conserving terms 
anymore. It can be written in terms of two structure functions, 

W^Bl + W'lon-A + W'lon-B = ^aZ^'S'J^ - XbZ^'S'Jl , (4-33) 



with 



2,1 



d^kTTz^giAxB, kl) + 2MbXb9t{xb) 



2M, 



b 



+ 



d^kTT^hliixA, kl) + 2MAXAhl{xA) 



Ma 

+2mg'l{xA)K{xB) 



hUxB) 



(4.34) 
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where the sum runs over quark and antiquark flavors. The other structure function, U2 1 , 
can be obtained from Eq. ( ^4.34| ) by replacing A ^ B. It is easily shown that they cannot 
be written in terms of grix) and hL{x) only. However, using Eqs. ( p.29|) and ( p.34[) , one can 
simultaneously eliminate the explicit fcr-integrals and quark-mass term; 



jjLT 



^ el[MB gl{xA)xB [grixB) + ~gT{xB 



+MAXA[hl{xA) + hl{xA)\hl{xB)]. 



(4.35) 



C. Discussion 



Having derived the DY hadron tensor, it is straightforward to calculate the cross section 
and double-spin asymmetries after contracting the leptonic tensor, written in terms of the 
angles 9 and 0, with the hadronic tensor in Eq. ( [4.33| ). The only nonzero 0{1/Q) contri- 
butions come from TL or LT scattering, of which we only consider the latter; the former 
simply follows by replacing A ^ B. The result is 



da{XA, Sbt) _ 
dxAdxBdQ 4(5^ 



'1 + cos^ e) Wt + Xa sin 29 cos( 



(4.36) 



neglecting O^l/Q"^) contributions. Here Wt = (1/3) /f (xyi)/^ (x^) is the leading- 
order unpolarized structure function. As transverse momentum has been integrated over, 
this result depends only on the relative azimuthal angle, 0—0^, between the lepton scattering 
plane and the transverse polarization vector. The asymmetry following from Eq. (^4.36|) reads 



A 



LT 



o'iXA, Sbt) 



*BT) 



cr{XA, Sbt) + o-{Xa, -Sbt) 



A, 



sin20cos(' 



^b) uZ 



1 + cos2 e 



(4.37) 



Taking = —1, 0b = 0, and writing out the structure functions, this becomes 



A 



sin 2^ cos 1 



LT 



1 + cos2 e Q 
j:ael{MBg'i{xA)xB 



9t{xb) + fxixB)] + MAXA[hl{xA) + hl{xA)]K{xB)} 



EaeimxA)mxB) 



(4.38) 



We emphasize again that Eq. ([4.38|) is the general result for the 0{1/Q) LT-asymmetry. 
No assumptions have been made for the fc-r-dependence of the distributions. If we assume 



5(fc|i) transverse momentum distributions and m = 0, the relations ( p.29| ) and (|2.34|) reduce 



to grix) = gT{x) and hiix) = hiix), respectively. Inserting these into Eq. ( [4.38| ), we arrive 
at the result derived by Jaffe and Ji in Ref. (taking Ma = Mb = M). 

The bag model gives an estimate for the quark distributions. Using the calculations 
of Refs. HJl^ (see Sec. H), we plot in Fig. ^ the (xa, a;B)-dependent part (without the 



angular dependence and the factor M/Q) of the asymmetry ( [4.381 ). It turns out that the 
zero-transverse-momentum zero-m approximation is of roughly the same form, but shifted 
upwards. This can also be seen from the cut along the line xa = xb, depicted in Fig. |12. 
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Clearly, the two results differ considerably. This could have been foreseen from the obser- 
vation, made in Sec. ||, that the zero-transverse-momentum zero-m relations grix) = grix) 
and = in the bag are severely violated. 

In sum, we have carefully analyzed the twist-three polarized distributions grix) and 
hiix) as to their dynamical content in a nonstandard (i.e., non-OPE) way. They are sen- 
sitive to nonzero intrinsic transverse momentum, as is also clear from the analysis of the 
physical processes; polarized deep-inelastic and Drell-Yan scattering at 0{1/Q). For the 
latter process we found that the LT-asymmetry, in addition to gT{,x) and hi^x), contains 
the quark-gluon-quark correlations gT^x) and hi^x). 
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FIGURES 



FIG. 1. The blob representing the quark-quark correlation function $q,^(PS'; k). 

FIG. 2. The quark-gluon-quark correlation function A4^^(PS'; k,p). 

FIG. 3. The physical area in the cir-plane for x = 1/2. The lower boundary is given by the 
line T = a — 1, the upper hy t = xa — x"^ . 

FIG. 4. Distribution function gxix) for a massless quark in the bag (solid line), its 
twist-two Wandzura-Wilczek part dy[gi{y)/y] (dotted) and its interaction-dependent part 
grix) — dy[gT{y)/y] (dashed) |5|. The dot-dashed line is grix) which diverges like l/x for 
X ^0. 

FIG. 5. h2 (solid curve) in the bag is twice g2 (dashed). 

FIG. 6. h^^x) for a massless quark in the bag (solid line), its twist-two part 2x dy[hi{y)/y'^] 
(dotted), and its interaction-dependent part hL{x) — 2x dy[hL{y) / y"^] (dashed) H^. The distri- 
bution hiix) (dot-dashed) diverges like l/x for x — > 0. 

FIG. 7. The quark Born diagram for DIS. 

FIG. 8. Gluon diagrams contributing to 0{l/Q) DIS. 

FIG. 9. DY quark Born diagram. 

FIG. 10. DY quark diagrams with one transverse-gluon insertion. 

FIG. 11. (x^, xb) -dependent part oi Alt [Eq. ( [4.38|) ] in the bag. 

FIG. 12. The bag-model (x^, XB)-dependent part of Alt for xa = xb (solid), and its ze- 
ro-transverse-momentum approximation (dashed) . 
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